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ABSTRACT 

A  full-spectral  third-generation  ocean  wind-wave  model,  Wavewatch-III,  has  been  implemented  in  the  South 
China  Sea  (SCS)  for  investigating  wind-wave  characteristics.  This  model  was  developed  at  the  Ocean  Modeling 
Branch  of  the  National  Centers  for  Environmental  Prediction  (NCEP).  The  NASA  QuickSCAT  data  (0.25° 
resolution)  2  times  daily  were  used  to  simulate  the  wind  waves  for  the  entire  year  of  2000.  The  significant  wave 
heights  from  Wavewatch-III  are  compared  to  the  TOPEX/Poseidon  (T/P)  significant  wave  height  data  over  the 
satellite  crossover  points  in  SCS.  The  model  errors  of  significant  wave  height  have  Gaussian-type  distribution 
with  a  small  mean  value  of  0.02  m  (almost  no  bias).  The  model  errors  are  comparable  to  the  T/P  altimeter 
accuracy  (0.5  m)  in  the  central  SCS  and  are  smaller  than  the  T/P  altimeter  accuracy  in  the  northern  and  southern 
SCS,  which  indicates  the  capability  of  Wavewatch-III  for  SCS  wave  simulation. 


1.  Introduction 

The  South  China  Sea  (SCS)  is  a  semienclosed  tropical 
sea  located  between  the  Asian  landmass  to  the  north 
and  west,  the  Philippine  Islands  to  the  east,  Borneo  to 
the  southeast,  and  Indonesia  to  the  south  (Fig.  1),  a  total 
area  of  3.5  X  106  km2.  It  connects  to  the  East  China 
Sea  (through  Taiwan  Strait),  the  Pacific  Ocean  (through 
Luzon  Strait),  the  Sulu  Sea,  the  Java  Sea  (through  Gas¬ 
per  and  Karimata  Straits),  and  to  the  Indian  Ocean 
(through  the  Strait  of  Malacca).  All  of  these  straits  are 
shallow  except  Luzon  Strait,  the  maximum  depth  of 
which  is  1800  m.  The  elliptically  shaped  central  deep 
basin  is  1900  km  along  its  major  axis  (northeast-south- 
west)  and  approximately  1100  km  along  its  minor  axis, 
and  extends  to  over  4000  m  deep. 

The  SCS  is  under  the  influence  of  monsoon  winds  and 
synoptic  systems  such  as  fronts  and  tropical  cyclones 
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(Chu  et  al.  1999a, b,  2000).  From  November  to  March, 
the  northeasterly  winter  monsoon  winds  correspond  to 
monthly  mean  January  2000  wind  speeds  of  near  10  m 
s  1  for  the  whole  SCS  (Fig.  2a).  From  April  to  August, 
the  southwesterly  summer  monsoon  winds  result  in 
monthly  mean  July  2000  wind  speeds  of  approximate  8 
m  s  1  in  the  southern  SCS  and  4  m  s  in  the  northern 
SCS  (Fig.  2b).  The  monthly  mean  winds  (Fig.  2)  are 
typical  for  monsoon  winds.  Highly  variable  winds  and 
surface  currents  are  observed  during  the  transitional  pe¬ 
riods.  Moreover,  synoptic  systems  often  pass  by  the  SCS 
and  cause  temporally  and  spatially  varying  wind  fields. 

A  fully  spectral  third-generation  ocean  wind-wave 
model,  Wavewatch-III  (henceforth  denoted  as 
WWATCH),  has  been  recently  developed  at  the  Ocean 
Modeling  Branch  of  the  Environmental  Modeling  Center 
of  the  National  Centers  for  Environmental  Prediction 
(NCEP)  for  the  regional  sea  wave  prediction.  It  was  built 
on  the  base  of  Wavewatch-I  and  Wavewatch-II  as  de¬ 
veloped  at  the  Delft  University  of  Technology  and  the 
National  Aeronautics  and  Space  Administration  (NASA) 
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Fig.  1.  Geography  and  isobaths  showing  the  bottom  topography  of  the  South  China  Sea. 


Goddard  Space  Flight  Center,  respectively  (Tolman 
1999). 

WWATCH  should  be  evaluated  before  the  practical 
use.  The  highly  variable  wind  systems  and  complicated 
topography  make  SCS  a  perfect  location  for  WWATCH 
evaluation.  Usually  the  in  situ  wind-wave  data  are  main¬ 
ly  collected  from  voluntary  ships  and  wave  buoys.  How¬ 
ever,  in  the  SCS  sparse  voluntary  ship  data  and  no  wave 
buoy  data  are  available.  Remote  sensing  is  an  important 
source  for  the  wind-wave  data.  Several  satellites  have 
been  launched  with  altimetry,  such  as  the  TOPEX/Po- 
seidon  (T/P)  satellite  and  the  European  Remote  Sensing 
Satellites  ERS-1/2.  The  subsequent  sections  describe  the 
WWATCH  evaluation  using  the  T/P  significant  wave 
height  (SWH)  data. 

The  outline  of  this  paper  is  as  follows.  A  description 
of  data  and  WWATCH  is  given  in  sections  2  and  3. 
Methodology  and  model  evaluation  are  given  in  sections 
4  and  5.  In  section  6,  we  present  our  conclusions. 

2.  Data 

a.  T/P  SWH  and  wind  speed 

The  T/P  satellite,  jointly  launched  by  NASA  and  the 
French  Space  Agency,  the  Center  National  d’Etudes  Spa- 
tiales  (CNES),  in  August  1992,  carried  a  state-of-the-art 


radar  altimetry  system  (Fu  et  al.  1994).  In  addition  to 
precise  measurements  of  the  distance  between  the  sat¬ 
ellite  and  the  surface,  SWH  and  wind  speed  are  derived 
from  the  shape  of  the  leading  edge  of  the  returning  radar 
pulse.  The  accuracy  of  SWH  measurement  by  T/P  was 
within  the  accuracy  of  the  Geosat  measurements  (Cal¬ 
lahan  et  al.  1994),  that  is,  10%  or  0.5  m,  whichever  is 
greater  (Dobson  et  al.  1987).  The  T/P  satellite  was  ma¬ 
neuvered  into  a  9.9156-day  repeat  period  during  which 
two  T/P  SWH  and  wind  speed  data  are  available  at  each 
crossover  point.  Time  series  of  SWH  at  20  crossover 
points  (Fig.  3a)  and  at  four  tracks  (051,  229,  152,  190 
in  Fig.  3b)  for  2000  are  used  to  evaluate  WWATCH. 

b.  QuikSCAT  sea  surface  winds 

NASA  launched  the  microwave  scatterometer 
SeaWinds  on  the  QuikBird  satellite  in  June  1999.  This 
instrument  is  referred  to  as  QuikSCAT  (or  QSCAT). 
QSCAT  is  essentially  a  radar  device  that  transmits  radar 
pulses  down  to  the  earth’s  surface  and  then  measures 
the  power  that  is  scattered  back  to  the  instrument.  This 
“backscattered”  power  is  a  measure  of  surface  rough¬ 
ness.  For  water  surfaces,  the  surface  roughness  is  highly 
correlated  with  the  near-surface  wind  speed  and  direc¬ 
tion.  Hence,  wind  speed  and  direction  at  a  height  of 
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Fig.  2.  Monthly  mean  wind  speed  at  10-m  height  computed  from  the  QSCAT  data:  (a)  Jan  and  (b)  Jul  2000. 


Fig.  3.  The  T/P  satellite  (a)  crossover  points  and  (b)  tracks  in  the  SCS. 
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Table  1.  Model  setting  for  this  study. 


Switch  parameter 


Characteristics 


DUM 

LRB8 

SHRD 

SEED 

GRD1 

SP1 

PR2 

ST2 

STAB2 

NL1 

BT1 

WIND2 

CUR2 

ol 

o2 

o2a 

o3 

o4 

o5 

06 

o7 


Dummy  to  be  used  if  WWATCH  is  to  be  installed  on  previously  untried  hardware 
8-byte  words 

Shared  memory  model,  no  message  passing 
Seeding  of  high-frequency  energy 

Settings  directly  hardwired  to  user-defined  spatial  grids  (spherical  coordinate  with  0.25°  grids) 
User-defined  spectral  grids 

Ultimate  quickest  propagation  scheme  with  Booij  and  Holthuijsen  (1987)  dispersion  correction 
Tolman  and  Chalikov  (1996)  source  term  package 

Enable  stability  correction  for  Tolman  and  Chalikov  (1996)  source  term  package 
Nonlinear  interaction  (DIA) 

JONS  WAP  bottom  friction  formulation 

Approximately  quadratic  interpolation 

Approximately  quadratic  interpolation 

Output  of  boundary  points  in  grid  preprocessor 

Output  of  the  gridpoint  status  map  in  grid  preprocessor 

Generation  of  landsea  mask  file  mask.ww3  in  grid  preprocessor 

Additional  output  in  loop  over  fields  in  field  preprocessor 

Print  plot  of  normalized  ID  energy  spectrum  in  initial  conditions  program 

2D  energy  spectrum 

Spatial  distribution  of  wave  heights  (not  adapted  for  distribution  memory) 

Echo  input  data  for  homogeneous  fields  in  generic  shell 


10  m  over  the  ocean  surface  are  retrieved  from  mea¬ 
surements  of  the  QSCAT  backscattered  power.  The 
backscattered  power  is  also  affected  by  rain. 

The  QSCAT  Level-3  global  surface  wind  dataset  con¬ 
sists  of  twice-daily  gridded  values  (0.25°  X  0.25°)  of 
scalar  wind  speed,  meridional  and  zonal  components  of 
wind  velocity,  wind  speed  squared,  and  time  given  in 
fraction  of  a  day.  The  presence  of  rain  may  degrade  the 
accuracy  of  the  QSCAT  winds.  Such  degradation  is  in¬ 
dicated  by  rain  probability  determined  using  the  multi¬ 
dimensional  histogram  (MUDH)  rain  flagging  technique. 
Data  are  currently  available  in  hierarchical  data  format 
(HDF)  and  exist  from  19  July  1999  to  the  present. 

3.  WWATCH  model 

a.  Description 

The  wave  spectrum  F  is  generally  a  function  of  all 
phase  parameters  (i.e.,  wavenumber  k,  direction  9,  in¬ 
trinsic  frequency  a,  and  absolute  frequency  a>),  space 
(x),  and  time  ( t ): 

F  —  F(k,  9,  cr,  co\  x,  t). 

However,  the  individual  spectral  components  are  usually 
assumed  to  satisfy  the  linear  wave  theory  (locally)  and 
to  follow  the  dispersion  relation, 

cr2  =  gk  tanhM  (1) 

co  =  a  +  k  •  U,  (2) 

where  d  is  the  mean  water  depth  and  U  is  the  (depth  and 
time  averaged)  current  velocity.  When  the  current  ve¬ 
locity  vanishes,  only  two-phase  parameters  among  (cr,  k, 
9)  are  independent.  Current  wave  models  use  the  fre¬ 
quency  direction  (cr,  9)  as  the  independent  phase  vari¬ 
ables. 


WWATCH  uses  the  wavenumber  direction  (k,  9)  as 
the  independent  phase  variables.  Without  currents,  the 
energy  of  a  wave  package  is  conserved.  With  currents 
the  energy  of  a  spectral  component  is  no  longer  conserved 
(Longuet-Higgins  and  Stewart  1961),  but  the  wave  action 
spectrum,  N(k,  9\  x,  t)  =  F(k,  9;  x,  f)/cr,  is  conserved 
(Whitham  1965;  Bretherthon  and  Garrett  1968).  In 
WWATCH,  the  balance  equation  is  for  the  wave  action 
spectrum. 

b.  Model  setting 

WWATCH  has  two  types  (mandatory  and  optional) 
of  model  switches  for  users  to  choose.  Table  1  lists  the 
model  setting  and  optional  switches  for  this  study.  For 
example,  spatial  and  spectral  grids  are  user  defined;  the 
ultimate  quickest  propagation  scheme  is  selected  with 
the  dispersion  correction  from  Booij  and  Holthuijsen 
(1987);  nonlinear  interactions  are  included;  and  the 
source  term  parameterization,  following  Tolman  and 
Chalikov  (1996),  consists  of  four  parts:  wind  input,  non¬ 
linear  wave-wave  interaction,  dissipation,  and  wave- 
bottom  interaction.  The  output  of  WWATCH  consists 
of  the  traditional  frequency-direction  spectrum  F{cr,  9), 
which  is  calculated  from  F(k,  9)  using  Jacobean  trans¬ 
formations. 

c.  Discretization 

The  model  is  implemented  for  SCS  (0°-25°N,  105°- 
122°E)  using  realistic  bathymetry  data  from  the  Naval 
Oceanographic  Office  DBDB5  database  and  a  regularly 
spaced  latitude-longitude  grid  with  the  grid  spacing 
0.25°  (i.e.,  AA  =  A</>  =  0.25°). 

The  wavenumber  grid  spacing  is  determined  by  the 
frequency  intervals  (total  25) 
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i  =  X„(Tm,  m  =  0,  1,  ,  24,  (3) 

with 

2Q  =  1.1,  o-0  =  0.0418.  (4) 

The  wave  direction  (6)  grid  spacing  is  15°  (i.e.,  AO  = 
15°). 

Four  time  steps  are  used  in  WWATCH  to  reach  com¬ 
putational  efficiency:  (a)  a  global  time  step  (300  s)  for  the 
propagation  of  the  entire  solution,  (b)  a  spatial  time  step 
(300  s)  representing  the  spatial  propagation,  (c)  a  spectral 
time  step  (300  s)  for  intraspectral  propagation,  and  (d)  a 
source  time  step  (100  s)  for  the  source  term  integration. 

d.  Wind  input  and  friction  velocities 

The  surface  winds  (W)  at  each  source  time  step 
(100  s)  for  the  year  2000  is  interpolated  from  NASA 
SeaWinds  on  twice-daily  QuikSCAT  (QSCAT)  Level- 
3  gridded  ocean  wind  vectors  with  0.25°  horizontal  res¬ 
olutions.  The  friction  velocities  are  needed  for  the  input 
source  function  Sin.  In  WWATCH,  the  friction  velocity 
(uj  is  computed  from  the  wind  speed  (W)  at  a  given 
reference  height  zr,  in  terms  of  a  drag  coefficient  C, 
(Tolman  and  Chalikov  1996): 

ul  =  C*W2(zr).  (5) 

An  iteration  process  is  adopted  to  obtain  uif.  Wu’s  (1982) 
empirical  relation 

CD(W)  =  1.2875  X  10-3  W  <  7.5  m  s”1, 

CD(W)  =  (0.8  +  0.065 W)  X  103  W  >  7.5  m  s"1, 

(6) 

is  used  to  obtain 


uf  =  VC ~W(z10),  (7) 

as  the  first-guess  friction  velocity.  Here,  z10  =  10  m. 
The  iteration  stops  when  the  change  of  the  friction  ve¬ 
locity  is  smaller  than  a  prescribed  criterion.  Such  iter¬ 
ations  are  performed  during  the  model  initialization  but 
are  not  necessary  during  the  actual  model  run,  as  u,,. 
changes  slowly  (Tolman  1999).  The  effect  of  the  at¬ 
mospheric  instability  on  the  friction  velocity  is  param¬ 
eterized  using  an  effective  wind  speed  We  (Tolman  and 
Booij  1998),  which  depends  on  the  surface  air  and  sea 
temperature  difference. 

e.  Model  integration 

WWATCH  is  integrated  with  twice-daily  gridded 
QSCAT  ocean  wind  vectors  (0.25°)  from  the  Joint  North 
Sea  Wave  Project  (JONSWAP)  1973  wave  spectra  (Has- 
selmann  et  al.  1980)  on  3  January  (no  sufficient  wind 
data  on  1-2  January  for  SCS)  until  31  December  2000. 
The  model  S  WH  data  are  interpolated  into  the  T/P  cross¬ 
over  points.  At  each  crossover  point,  there  are  M  pairs 
(approximately  72)  of  modeled  (Hm)  and  observed  (HJ 
SWH  data  in  2000  (around  two  pairs  per  10  days). 

4.  Methodology  of  verification 

a.  Verification  at  crossover  points 

The  difference  of  the  modeled  and  observed  SWH, 

A H  =  HJx,  y,  t )  -  H fix,  y,  t),  (8) 

represents  the  model  error.  Bias,  root-mean-square  error 
(rmse),  and  correlation  coefficient  (cc)  for  each  cross¬ 
over  point 


bias(x,  y) 


rmse(x,  y) 


cc(x,  y) 


1  M 

M  S  *  ''>• 

fl  " 

J  g  [M/C.  ,,  .,)]■, 

M 

X  [HJx,  y,  t,)  -  HJx,  y)][Hfix,  y,  t,)  -  Hfix,  y)] 


/2  y,  t,)  -  HJx,  y)]2  /2  [Hfix,  y,  f,)  -  Hfix,  y)]2 


(9) 

(10) 


(11) 


are  used  to  verify  WWATCH.  Here  HJx,  y)  and 
Hfix,  y)  are  temporal  mean  modeled  and  observed  SWH, 

|  M 

HJx,  y)  =  —  2  HJx,  y >  tfi 

M  i=i 

j  M 

H„(x,  y)  =  —  2  Hfix,  y,  tfi 

lVl  i=  l 


at  the  crossover  points.  The  T  value, 

cc  Vm  -  2 
T  Vl  -  cc2  ’ 


(13) 


(22)  with  the  degrees  of  freedom  of  (M  —  2)  is  used  for  the 
significant  test  of  the  correlation  coefficient. 
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Fig.  4.  Simulated  monthly  mean  SWH  using  WWATCH:  (a)  Jan  and  (b)  Jul  2000. 


b.  Verification  at  time  instance 
Bias  and  rmse  for  time  instance  t 

bias(r)  =  ^  2  yk,  t),  (14) 

i,k 


rmse(f)  =  hz  2  yk,  t)]2,  (15) 

yiy  j.k 

are  also  used  to  verify  WWATCH. 

5.  Model  results 

Two  types  of  model-observation  comparison  are  con¬ 
ducted:  1)  horizontal  SWH  pattern  using  monthly  mean 
data  and  2)  model  error  statistics  using  synoptic  crossover 
point  SWH  data  (model  and  observation).  To  show  the 
value  added  by  using  WWATCH,  the  Pierson-Moskowitz 
(P-M)  spectrum  (Pierson  and  Moskowitz  1964)  is  also 
used  for  comparison  in  monthly  mean  horizontal  SWH 
pattern. 

a.  Monthly  mean  SWH 

Three  sets  of  monthly  mean  data  are  calculated  from 
(a)  modeled  SWH  using  WWATCH,  (b)  calculated  SWH 
using  the  P-M  spectrum  with  the  same  QSCAT  winds, 
and  (c)  observed  SWH  from  T/P.  A  common  feature  in 
the  simulated  (Fig.  4),  calculated  (Fig.  5),  and  observed 
(Fig.  6)  data  shows  a  higher  SWH  in  January  2000  than 
in  July  2000. 

In  January  2000,  a  southwest-to-northeast-oriented 
high-SWH  region  (>2.0  m)  is  comparable  (north  of  5°N) 


in  the  WWATCH  simulation  (Fig.  4a)  and  in  the  T/P  data 
(Fig.  6a).  However,  this  high-SWH  region  is  split  into 
two  smaller  ones  in  the  calculated  (from  P-M  spectrum) 
field,  with  a  major  one  occurring  north  of  15°N  and  a 
minor  one  near  the  southern  Vietnamese  coast  (Fig.  5a). 
The  area  with  SWH  larger  than  2.5  m  in  the  WWATCH 
simulation  (15°-20°N,  1 13°-1 17°E)  is  comparable  to  that 
in  the  T/P  data  (13°-20°N,  112°-117°E).  Due  to  its  in¬ 
dependence  on  fetch,  the  P-M  results  show  the  SWH 
maximum  of  3  m  near  the  northeast  boundaries  (upwind), 
which  is  not  found  in  the  simulated  (WWATCH)  and 
observed  fields. 

In  July  2000,  the  mean  SWHs  are  higher  in  the  north¬ 
ern  and  central  SCS  (north  of  10°N)  than  in  the  southern 
SCS  (south  of  10°N),  with  values  larger  than  2.25  m  in 
the  WWATCH  simulated  field  (Fig.  4b)  and  larger  than 
2.5  m  in  the  calculated  (from  the  P-M  spectrum)  (Fig. 
5b)  and  observed  fields  (Fig.  6b).  However,  the  maximum 
SWH  values  are  located  at  1 1°-15°N,  1 15°-120°E  in  sim¬ 
ulated  and  calculated  fields  and  at  15°-20°N,  1 13°-1 16°E 
in  the  observed  field. 

Due  to  its  independence  on  fetch,  the  P-M  results  show 
the  SWH  maximum  of  3  m  near  the  northeast  boundaries 
(upwind)  in  January  and  of  2  m  near  the  southern  bound¬ 
aries  (upwind)  in  July,  which  is  not  found  in  the  simulated 
(WWATCH)  and  observed  fields.  WWATCH  simulates 
the  seasonal  variability  of  SWH  reasonably  well.  SWH 
is  larger  in  the  winter  than  in  the  summer  monsoon  sea¬ 
son.  The  orientation  of  the  high-SWH  region  coincides 
with  the  orientation  of  the  monsoon  winds  (Fig.  2).  The 
observed  data  show  some  more  intensification  in  the 
maximum  values.  The  SWHs  show  more  irregular  pat- 
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(a)  P-M  SWH(m)  in  JAN  <b)  P-M  SWH(m)  in  JUL 


Fig.  6.  Monthly  mean  SWH  from  T/P  data:  (a)  Jan  and  (b)  Jul  2000. 
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(a) 


AH(m) 


Fig.  7.  Model  accuracy  statistics:  (a)  histogram  of  model  error  and 
(b)  scatter  diagram  of  modeled  (Hm)  and  observed  (Ha)  SWH  for  all 
the  crossover  points. 


terns  using  the  P-M  spectrum  than  using  WWATCH. 
This  is  due  to  the  coarse  resolution  in  time  of  the  wind 
input  (twice  daily)  that  smoothes  out  extremes  such  as 
storm  event  (Chu  et  al.  2000;  Chen  et  al.  2002). 

b.  Statistical  evaluation 

The  model-observation  comparison  is  conducted  at  all 
20  crossover  points  in  the  SCS.  Each  one  contains  ap¬ 
proximate  72  pairs  of  modeled  ( Hm )  and  observed  (Ha) 
SWH  data  in  2000.  The  total  number  of  pairs  is  1330. 
The  difference  between  the  two,  A H  =  Hm  —  Ha,  is 


considered  as  the  model  error.  The  histogram  of  AH  = 
[=  Hm  —  //„]  (Fig.  7a)  for  all  20  crossover  points  shows 
a  Gaussian-type  distribution  with  mean  value  (—0.02  m) 
and  with  comparable  sample  numbers  of  positive  AH 
(633)  with  negative  AH  (697).  The  scatter  diagrams  for 
Hm  and  Ha  show  clustering  of  points  approximately 
around  the  line  of  Hm  =  Ha  (Fig.  7b).  The  rmse  and  cc 
between  Hm  and  If,,  are  0.48  m  and  0.90,  respectively. 

The  scatter  diagrams  for  Hm  and  Ha  at  each  crossover 
point  (Fig.  8)  show  spatial  variability  of  the  error  statis¬ 
tics.  The  rmse  increases  from  0.34  m  at  the  southwest 
corner  of  the  SCS  near  Natuna  Island  [point  1  (2.01°N, 
106.3 1°E)]  to  0.95  m  in  the  central  SCS  near  Zhongsha 
Island  [point  15  (17.18°N,  114.81°E)]  and  west  of  Luzon 
[point  16  (17.20°N,  117.65°E)].  The  bias  varies  from 
-0.45  m  in  the  Gulf  of  Tonkin  [point  17  (20.59°N, 
107.73°E)]  to  0.33  m  at  the  southern  SCS  near  Nansha 
Island  [point  7  (9.8°N,  114.81°E)].  The  correlation  co¬ 
efficient  increases  from  0.55  in  the  Gulf  of  Tonkin  (point 
17)  to  0.95  in  the  central  SCS  near  Zhongsha  Island  (point 
15). 

Contours  of  bias,  rmse,  and  cc  for  the  entire  year  of 
2000  are  plotted  (Fig.  9)  to  understand  the  spatial  error 
variability.  A  positive  bias  occupies  a  large  portion  of 
the  SCS.  The  zero-bias  contour  follows  200-m  ba¬ 
thymetry  (Fig.  1)  with  negative  bias  on  the  continental 
shelf  (west  of  the  zero-bias  contour)  and  positive  bias 
in  the  deep  basin  (east  of  the  zero-bias  contour).  A 
negative  bias  larger  than  —0.4  m  is  found  in  the  Gulf 
of  Tonkin  and  a  positive  bias  larger  than  0.3  m  is  lo¬ 
cated  near  Nansha  Island  (10°N,  115°E)  (Fig.  9a).  This 
indicates  that  WWATCH  overpredicts  the  SWH  slight¬ 
ly  except  on  the  shallow  continental  shelf.  The  rmse 
of  SWH  is  above  0.5  m  in  the  central  SCS,  with  a 
maximum  larger  than  0.6  m  west  of  Palawan  (Fig.  9b). 
The  value  of  rmse  decreases  from  the  central  SCS  to 
the  other  two  regions  and  is  smaller  than  0.5  m  in  most 
of  the  northern  (west  of  Luzon)  and  southern  (south 
of  5°N)  SCS. 

Strong  negative  bias  in  the  Gulf  of  Tonkin  is  ob¬ 
tained  on  the  basis  of  observations  at  a  single  crossover 
point  (i.e.,  point  17),  which  is  quite  close  to  the  coast. 
The  T/P  altimeter  data  are  not  as  valuable  as  in  the 
deep  sea.  Recently,  two  high-kinetic -energy  zones 
were  identified  as  tidal  fronts  in  the  Qiongzhou  Strait 
and  near  the  southwestern  coast  of  the  Hainan  Island 
(Hu  et  al.  2003).  Tidal  effect  is  not  available  in  the 
model  simulation. 

The  cc  of  SWH  (Fig.  9c)  between  modeled  and 
T/P  data  in  2000  is  larger  than  0.85  almost  everywhere 
in  the  SCS,  except  in  the  Gulf  of  Tonkin.  The  T  value 
computed  using  (13)  for  cc  =  0.85,  M  =  72  is  T  = 
13.50.  For  confidence  coefficient  (1  —  a)  =  0.095,  the 
t  distribution  for  the  degree  of  freedom  of  statistics  for 
(M  —  2  =  70)  is  2.756  >  f0005  >  2.576.  Since  T(= 
13.50)  is  larger  than  t00 05,  the  correlation  coefficient 
between  modeled  and  T/P  SWH  data  is  significant. 
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Point  17:  107.7^E  20.59°N  Point  18:  113.4?E  20.58“N 


Point  19:  116.2^E  20.62°N  Point  20:  119.07°E  20,58°N 


Fig.  8.  Scatter  diagrams  of  SWHs  at  crossover  points  in  the  SCS. 


c.  Spatial  error  variability  along  three  ascending 
tracks 

Understanding  the  spatial  error  variability  during  dif¬ 
ferent  periods,  four  tracks  (051,  229,  153,  190;  see  Fig. 
3b)  covering  the  SCS  and  three  cycles  (270,  288,  302) 


are  selected  for  the  model  evaluation.  For  each  cycle, 
the  T/P  SWHs  and  wind  speeds  along  the  tracks  are 
compared  to  the  model-simulated  SWHs  and  QSCAT 
wind  speeds  (interpolated  with  the  same  temporal  res¬ 
olution  with  the  T/P  SWH  and  wind  speed  data). 
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(a)  (b)  (c) 


Fig.  9.  Distributions  of  SWH  (a)  bias,  (b)  rms  error,  and  (c)  correlation  coefficient  between  WWATCH  and  T/P  altimeter  data  for  the 

entire  year  of  2000. 


1)  Winter  monsoon  season  (cycle  270) 

During  cycle  270,  the  T/P  satellite  passed  over  tracks 
051,  153,  190,  and  229  on  14,  18,  19,  and  21  January, 
respectively.  The  simulated  and  observed  SWHs  on  the 
four  tracks  are  presented  in  the  left-hand  panels  of  Fig. 
10.  The  QSCAT  and  T/P  wind  speeds  on  the  four  tracks 
are  presented  in  the  right-hand  panels. 

Along  tracks  051,  229,  and  190  (eastern  and  middle 
SCS)  the  simulated  SWH  coincides  with  T/P  SWH  quite 
well  south  of  12°N,  with  errors  usually  less  than  0.3  m, 
and  deviates  from  T/P  SWH  north  of  12°N,  with  errors 
larger  than  1.0  m.  Along  track  153  (western  SCS)  the 
simulated  SWH  is  quite  different  from  T/P  SWH,  with 
errors  larger  than  1.0  m  except  for  the  10°-18°N  seg¬ 
ment. 

After  the  gridded  QSCAT  winds  are  interpolated  tem¬ 
porally  at  each  source  time  step  (100  s),  the  QSCAT 
and  T/P  wind  speed  pairs  can  be  easily  obtained  along 
the  T/P  tracks  with  no  spatial  and  temporal  lags.  QSCAT 
wind  speeds  agree  with  T/P  wind  speeds  quite  well 
along  the  tracks  where  the  wind  speed  does  not  have 
strong  spatial  variability,  such  as  along  track  190. 
QSCAT  wind  speeds  have  a  large  discrepancy  along 
track  153  south  of  10°N  and  track  229  north  of  20°N, 
where  the  wind  speed  has  strong  spatial  variability,  with 
a  difference  up  to  around  3-4  m  s  1 .  Comparison  be¬ 
tween  left-  and  right-hand  panels  shows  that  discrepancy 
in  SWH  is  correlated  with  wind  speed  error,  except 
along  track  229,  where  the  model  errors  do  not  well 
correlate  to  the  wind  errors.  This  may  imply  the  im¬ 
portance  of  other  source  functions  such  as  dissipation 
and  nonlinear  effect. 


2)  Summer  monsoon  season  (cycle  288) 

During  cycle  288,  the  T/P  satellite  passed  over  tracks 
051,  153,  190,  and  229  on  10,  14,  16,  and  17  July, 
respectively.  The  simulated  and  observed  SWHs  on  the 
four  tracks  are  presented  in  the  left-hand  panels  of  Fig. 
1 1 .  The  QSCAT  and  T/P  wind  speeds  on  the  four  tracks 
are  presented  in  the  right-hand  panels.  Along  tracks  153 
and  229  the  simulated  SWH  coincides  with  T/P  SWH 
quite  well,  with  errors  less  than  0.5  m  in  most  of  lati¬ 
tudes.  Along  tracks  051  and  190  the  simulated  SWH  is 
quite  different  from  T/P  SWH,  with  errors  larger  than 
1.0  m  in  most  latitudes. 

QSCAT  wind  speeds  agree  with  T/P  wind  speeds 
quite  well  along  track  153,  which  leads  to  a  better  SWH 
simulation  along  that  track.  QSCAT  wind  speeds  have 
a  large  discrepancy  from  the  T/P  wind  speeds  along 
track  153,  which  leads  to  larger  SWH  errors  along  that 
track.  Thus,  the  summer  monsoon  winds  seem  more 
difficult  to  resolve  in  the  model.  Quality  of  the  wind 
input  is  important  for  SWH  prediction. 

3)  Tropical  storm  rumbia  2000  (cycles  302 
and  303) 

Cycles  302  and  303,  during  which  Tropical  Storm 
Rumbia  2000  was  passing  SCS  with  maximum  guess 
strength  of  25-30  m  s_1,  are  selected  for  evaluation  of 
model  capability  under  strong  and  highly  varying  wind 
forcing.  Tropical  Storm  Rumbia  initially  formed  about 
800  km  east  of  the  Philippine  island  of  Mindanao  on 
28  November  2000  (Fig.  12).  After  formation,  Rumbia 
slowly  intensified  as  it  tracked  westward  through  the 
Philippine  Sea  toward  the  central  Philippine  Islands.  On 
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Cycle  270  Pass  190  19-Jan-2000  20:33:33(UTC) 


Cycle  270  Pass  229  21 -Jan-2000  09:14:28(UTC) 
20 


Fig. 


10  15 

Latitude 

10.  Comparison  between  (left)  modeled  (dotted)  and  T/P  observed  (circle)  SWHs  and  (right)  QSCAT  (dotted)  and 
T/P  observed  wind  speeds  (circle)  along  the  four  tracks  during  cycle  270  (14-21  Jan  2000). 


1  December,  it  made  landfall  over  the  eastern  Philip¬ 
pines  and  intensified  to  tropical  storm  strength,  18 
m  s_1  (measured  by  QSCAT).  It  continued  moving  west¬ 
erly  through  the  Philippine  Islands,  intensifying  at  a 
slow  rate.  After  it  entered  the  SCS,  it  continued  to  move 
westward,  and  weakened  and  disappeared  on  9  Decem¬ 
ber  2000  near  the  southern  tip  of  Vietnam. 


The  T/P  satellite  passed  over  track  190  (2  December) 
and  track  229  (3  December)  during  cycle  302,  and  over 
track  051  (6  December)  and  track  153  (10  December) 
during  cycle  303.  The  simulated  and  observed  SWHs 
(from  T/P  altimeter)  on  the  four  tracks  are  presented  in 
the  left-hand  panels  of  Fig.  13.  The  QSCAT  and  T/P 
wind  speeds  on  the  four  tracks  are  presented  in  the  right- 
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Cycle  288  Pass  051  10- Jul-2000  22:01:25(UTC) 


Cycle  288  Pass  153  14- Jul-2000  21:35:01(UTC) 


Cycle  288  Pass  190  16- Jul-2000  08:07:24(UTC) 


Cycle  288  Pass  229  17-Jul-2000  20:47:52(UTC) 


Fig.  11.  Comparison  between  (left)  modeled  (dotted)  and  T/P  observed  (circle)  SWHs  and  (right)  QSCAT  (dotted)  and 
T/P  observed  wind  speeds  (circle)  along  the  four  tracks  during  cycle  288  (10-17  Jul  2000). 


hand  panels.  Along  track  190  the  simulated  SWH  co¬ 
incides  with  T/P  SWH  quite  well  south  of  17°N,  with 
mean  error  less  than  0.5  m.  Along  track  229  the  sim¬ 
ulated  SWH  is  quite  different  from  T/P  SWH  in  the 
central  SCS  (5°-18°N).  QSCAT  wind  speeds  agree  with 
T/P  wind  speeds  quite  well  along  track  190  south  of 


18°N,  which  leads  to  a  better  SWH  simulation  along 
that  track.  Wind  speeds  are  quite  different  between 
QSCAT  and  T/P  in  the  central  SCS  (8°-18°N)  along 
track  229,  which  leads  to  larger  SWH  simulation  errors 
along  track  229. 

Poor  performance  during  Tropical  Storm  Rumbia  may 
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Tropical  Storm  33W  (Rumbia)  Best  Track  27  Nov  -  09  Dec  2000 


110°E  120°E  130°E 


Fig.  12.  Track  of  Tropical  Storm  Rumbia  (2000). 

be  caused  by  too-coarse  grid  resolution  and  lack  of  at¬ 
mosphere-wave-ocean  coupling.  Chen  et  al.  (2002) 
pointed  out  that  simulation  of  typhoon  or  hurricane 
needs  resolution  from  1/6°  for  general  structure  to  2  km 
for  detail  eyewall  structure.  With  the  extreme  high 
winds,  intense  rainfall,  large  ocean  waves,  and  copious 
sea  spray  the  surface  flux  parameterization  should  be 
updated.  This  is  illustrated  by  numerical  simulation  of 
air-sea  interaction  under  high  wind  conditions  using  a 
coupled  atmosphere-ocean-wave  modeling  system 
(Bao  et  al.  2000).  Improvement  of  WWATCH  for  high 
wind  conditions  needs  finer  resolution  and  more  realistic 
parameterization  for  surface  momentum  flux. 

d.  Temporal  error  variability  at  selected  crossover 

points 

Understanding  the  temporal  error  variability,  time  se¬ 
ries  of  Hm  and  //„  at  four  crossover  points  are  presented 
(Fig.  14).  WWATCH  shows  the  capability  in  simulating 
SWH.  Crossover  point  2  (Fig.  14a)  represents  the  south¬ 
ern  SCS  (Fig.  3).  Bias,  rmse,  and  cc  at  that  point  are 
0.12  m,  0.41  m,  and  0.92,  respectively.  The  model  errors 
do  not  have  evident  seasonal  variability;  that  is,  the  model 
performance  is  equivalent  during  the  prevailing  monsoon 
seasons  and  during  the  monsoon  transition  periods. 

Crossover  point  6  represents  the  central  SCS  near  10°N 
(Fig.  3).  Bias,  rmse,  and  cc  at  that  point  are  0.03  m,  0.55 
m,  and  0.90,  respectively.  The  model  errors  (Fig.  14b) 
are  smaller  during  the  (northeast)  winter  monsoon  season 
(November-March)  than  during  the  (southwest)  summer 
monsoon  season  (May-September). 

Crossover  point  1 1  represents  the  eastern  part  of  the 
central  SCS  (15°N)  near  Luzon  Island  (Fig.  3).  Bias, 
rmse,  and  cc  at  that  point  are  0.09  m,  0.49  m,  and  0.90, 
respectively.  Different  from  crossover  point  6,  the  model 
errors  (Fig.  14c)  are  larger  during  the  (northeast)  winter 
monsoon  season  (November-March)  than  during  the 
(southwest)  summer  monsoon  season  (May-September). 
For  example,  the  SWHs  are  all  well  predicted  during  the 


Table  2.  Root-mean-square  difference  and  maximum  difference 
between  QSCAT  and  T/P  wind  speeds. 


Cycle 

No. 

Track  No. 

No.  of  data 
pairs 

Rms 

difference 
(m  s-1) 

Max  difference 
(m  s-1) 

270 

051 

372 

1.12 

2.87 

270 

153 

405 

2.19 

4.47 

270 

190 

353 

2.42 

4.43 

270 

229 

415 

1.90 

4.46 

288 

051 

333 

2.83 

9.98 

288 

153 

401 

1.49 

4.32 

288 

190 

341 

2.40 

8.00 

288 

229 

404 

2.84 

8.73 

302 

190 

332 

1.52 

11.28 

302 

229 

426 

3.70 

11.92 

303 

051 

370 

2.30 

6.72 

303 

153 

380 

2.36 

6.79 

summer  monsoon  season;  however,  during  the  winter 
monsoon  period  on  2  February  the  model-simulated 
SWH  (3.27  m)  is  much  larger  than  the  observed  SWH 
(2.10  m). 

Crossover  point  19  represents  the  north  SCS  near  20°N 
(Fig.  3).  Bias,  rmse,  and  cc  at  that  point  are  0.07  m,  0.43 
m,  and  0.94,  respectively.  The  model  errors  (Fig.  14d) 
are  larger  during  the  (northeast)  winter  monsoon  season 
(November-March)  than  during  the  (southwest)  summer 
monsoon  season  (May-September).  For  example,  the 
SWHs  are  all  well  predicted  during  the  summer  monsoon 
season;  however,  during  the  winter  monsoon  period  on 
7  December  the  model-simulated  SWH  (4.61  m)  is  much 
larger  than  the  observed  SWH  (3.50  m). 

e.  Temporal  error  variability  for  the  whole  SCS 

The  monthly  mean  bias  and  rmse  averaged  over  all 
of  the  crossover  points  in  the  SCS  are  presented  in  Fig. 
15.  They  demonstrate  the  temporal  error  variability  for 
the  whole  SCS  tracks.  WWATCH  has  very  low  bias 
(—0.01  to  0.04  m)  in  predicting  SWH,  with  a  maximum 
(positive  bias)  value  of  0.04  m  in  March  and  a  minimum 
(negative  bias)  value  of  —0.01  m  in  April.  The  rmse 
has  a  minimum  value  of  0.39  m  in  March  and  a  max¬ 
imum  value  of  0.48  m  in  December. 

6.  Uncertainty  in  surface  wind  data 

Discrepancies  between  QSCAT  and  T/P  wind  speeds 
(right-hand  panels  of  Figs.  10-13)  are  found  especially 
during  the  Tropical  Storm  Rumbia  period  (1-9  Decem¬ 
ber  2000),  where  the  QSCAT  winds  are  stronger  than 
the  T/P  winds.  Table  2  lists  detail  comparison  between 
the  two  wind  datasets.  The  root-mean-square  difference 
between  the  two  wind  datasets  increases  from  a  mini¬ 
mum  value  of  1.12  m  s_1  on  14  January  2000  (cycle 
270,  track  051)  to  a  maximum  value  of  3.70  m  s_1  on 
3  December  2000  (cycle  302,  track  229).  The  maximum 
difference  between  the  two  wind  datasets  increases  from 
a  minimum  value  of  2.87  m  s_1  on  14  January  2000 
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Fig.  13.  Comparison  between  (left)  modeled  (dotted)  and  T/P  observed  (circle)  SWHs  and  (right)  QSCAT  (dotted)  and 
T/P  observed  wind  speeds  (circle)  along  the  four  tracks  during  cycles  302  and  303  (3-9  Dec  2000),  when  Tropical  Storm 
Rumbia  passed  over  the  SCS. 


(cycle  270,  track  051)  to  a  maximum  value  of  11.92 
m  s_1  on  3  December  2000  (cycle  302,  track  229).  Since 
the  QSCAT  winds  are  sensitive  to  rain,  the  large  dis¬ 
crepancy  during  the  Tropical  Storm  Rumbia  (2000)  pe¬ 
riod  may  be  caused  by  heavy  rain. 


7.  Conclusions 

Comparing  the  South  China  Sea  significant  wave 
height  hindcast  using  a  third-generation  wave  model 
(Wavewatch-III)  with  significant  wave  height  measured 
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Fig.  14.  Comparison  between  WWATCH  modeled  (solid  curve)  and  T/P  observed  (denoted  by  circle)  significant 
wave  heights  at  (a)  point  2.  (b)  point  6,  (c)  point  1 1,  and  (d)  point  19. 


by  TOPEX/Poseidon  altimeter  for  2000,  several  char¬ 
acteristics  of  the  model  errors  are  obtained  for  the  three 
subregions:  the  central,  northern,  and  southern  SCS. 

1)  Wavewatch-III  simulates  the  seasonal  variability  of 
SWH  reasonably  well  compared  to  the  T/P  SWH 
data.  July  2000  SWHs  are  higher  in  the  northern  and 
central  SCS  (north  of  10°N)  than  in  the  southern  SCS 
(south  of  10°N),  with  values  up  to  2.5  m  in  the 
WWATCH  simulated  field  and  3.0  m  in  the  observed 


fields.  The  orientation  of  the  high-SWH  region  co¬ 
incides  with  the  orientation  of  the  monsoon  winds. 

2)  The  model  errors  for  SWH  hindcast  have  Gaussian- 
type  distribution  with  mean  values  of  0.02  m  and 
with  slightly  more  sample  numbers  on  the  negative 
side  (697)  than  on  the  positive  side  (633).  The  root- 
mean-square  error  and  correlation  coefficient  be¬ 
tween  modeled  and  observed  significant  wave 
heights  are  0.48  m  and  0.90,  respectively. 

3)  The  model  errors  of  WWATCH  have  spatial  vari- 


rmse(m)  bias(m) 
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Fig.  15.  Temporal  evolution  of  (a)  bias  and  (b)  rmse  for  the 
whole  SCS. 


ability  with  overprediction  of  the  SWH,  except  on  the 
shallow  continental  shelf.  The  rmse  of  SWH  is  above 
0.5  m  in  the  central  SCS,  with  a  maximum  larger  than 
0.6  m  west  of  Palawan.  The  value  of  rmse  northward 
and  southward  decreases  from  the  central  SCS  and  is 
smaller  than  0.5  m  in  most  of  northern  (west  of  Luzon) 
and  southern  (south  of  5°N)  SCS. 

4)  Over  the  whole  SCS,  WWATCH  has  very  low  bias 
(—0.01  to  0.04  m)  in  predicting  SWH,  with  a  max¬ 
imum  (positive  bias)  value  of  0.04  m  in  March  and 
a  minimum  (negative  bias)  value  of  —0.01  m  in 
April.  The  root-mean-square  error  has  a  minimum 
value  of  0.39  m  in  March  and  a  maximum  value  of 
0.48  m  in  December. 

5)  The  model  errors  are  comparable  to  the  T/P  altimeter 
accuracy  (0.5  m)  in  the  central  SCS  and  smaller  than 
the  T/P  altimeter  accuracy  in  the  northern  and  south¬ 
ern  SCS,  which  indicates  the  capability  of  Wave- 
watch-III  for  SCS  wave  simulation. 

6)  Ocean  wave  models  are  mainly  forced  by  the  surface 
winds.  The  difference  between  QSCAT  and  T/P  wind 
speeds  in  this  study  indicates  evident  wind  uncer¬ 
tainty.  Modeling  errors  are  caused  by  uncertain  mod¬ 
el  external  forcing  (such  as  winds)  or  uncertain  mod¬ 
el  internal  structure  (such  as  imperfect  physics  and 
resolution).  Questions  arise.  What  is  the  contribution 
of  external  and  internal  uncertainty  on  modeling  er¬ 
rors?  How  does  error  propagate  from  winds  to  ocean 
waves?  Will  the  wind  error  be  amplified  or  damped 
after  it  enters  the  ocean  wave  models?  We  will  an¬ 
swer  these  questions  in  the  near  future. 
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